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Abstract. Ground-state and finite-temperature properties of the mixed spin-^ and 
spin-S 1 Ising-Heisenberg diamond chains are examined within an exact analytical 
approach based on the generalized decoration- iteration map. A particular emphasis is 
laid on the investigation of the effect of geometric frustration, which is generated by 
the competition between Heisenberg- and Ising-type exchange interactions. It is found 
that an interplay between the geometric frustration and quantum effects gives rise to 
several quantum ground states with entangled spin states in addition to some semi- 
classically ordered ones. Among the most interesting results to emerge from our study 
one could mention a rigorous evidence for quantized plateux in magnetization curves, 
an appearance of the round minimum in the thermal dependence of susceptibility 
times temperature data, double-peak zero-field specific heat curves, or an enhanced 
magnetocaloric effect when the frustration comes into play. The triple-peak specific 
heat curve is also detected when applying small external field to the system driven by 
the frustration into the disordered state. 
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1. Introduction 

Over the last three decades, the low- dimensional quantum spin models with competing 
(frustrated) interactions have attracted considerable research interest especially due 
to their extraordinary diverse ground-state behaviour. Geometrically frustrated spin 
systems constitute a special sub-class of the frustrated models that can be distinguished 
by incapability of spins, inherent in their lattice positions, to simultaneously minimize 
the ground-state energy of each individual spin-spin interaction pQ. As a rule, the 
quantum spin systems affected by a rather strong geometric frustration often exhibit an 
exotic non-magnetic ground state (which does not have its classical analogue) in addition 
to a rich variety of the semi-classically ordered ones [2]. It is worthy to notice, moreover, 
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that a subtle interplay between the geometric frustration and quantum fluctuations 
gives rise to a number of intriguing phenomena including the macroscopic degeneracy 
of ground state [Sj, order by disorder effect chirality [5 3 , quantum phase transitions 
[S|-[E], quantized plateaux in the magnetization curves 0-1101, double-peak specific heat 
curves [TTJ-[T31 , enhanced magnetocaloric effect [Hj-[IE], etc. 

Despite a significant amount of effort, there are only few frustrated spin-^ quantum 
Heisenberg models, such as Majumdar- Ghosh model on double chain J7|-[2n|, sawtooth 
(A) chain [21]- (21], or Shastry-Sutherland model [2H], for which precise analytic solution 
is available leastwise for the ground state. Nevertheless, it should be pointed out that 
frustrated quantum systems are in general rather difficult to deal with, since extensive 
numerical methods must be used in order to obtain a reliable estimate of their magnetic 
properties. From this point of view, the one- dimensional (ID) frustrated spin systems 
are the simplest systems with respect to accurate treatment. Of these systems, the spin- 
al quantum Heisenberg model with diamond chain topology is currently actively engaged 
in the investigation of geometric frustration. Interestingly, this simple quantum system 
turned out to have rather complex ground state; apart from the usual ferrimagnetic 
phase there are in fact several quantum dimerized and plaquette states involved in the 
zero-field ground-state phase diagram j2H]-[2El- Further studies devoted to the spin-| 
quantum Heisenberg model on the diamond chain have provided fairly accurate results 
for the ground-state phase diagram in a presence of the external field [20]- [SI], the spin 
gap 32j, the magnetization and susceptibility [33]. Another remarkable finding relates 
to the observation of an inversion phenomenon, which can be induced in the frustrated 
diamond chain through the exchange anisotropy [33JESj. Note that the ground state 
and thermodynamics of the mixed-spin diamond chains containing also higher-spin sites 
have already been particularly examined as well 36J- 40 . 

It is worthwhile to remark that ID frustrated spin systems have initially been 
introduced purely as toy models suitable for investigating the effect of spin frustration. 
However, a recent progress achieved in the design and controlled synthesis of molecular- 
based magnetic materials afforded another stimulus for testing ID frustrated spin 
systems by overcoming the lack of desirable model compounds. As a matter of fact, 
the rapidly expanding field of molecular engineering have led to fabrication of several 
coordination polymers, which can be regarded as genuive examples of the frustrated 
spin models (Hj. With the help of structural data known long ago [12], Kikuchi 
and co-workers [E]-[EI] have recognized Cu3(C03)2(OH) 2 (azurite) as an appropriate 
candidate for the diamond chain compound. Experimental data measured for the high- 
field magnetization, susceptibility, specific heat, NMR [I3j-[III| and ESR data |Tf]-|5T] 
have indeed confirmed Kikuchi's conjecture and it is now quite well established that the 
azurite represents the actual material for the frustrated diamond chain. It should be 
also mentioned that sufficiently strong frustration found in the azurite clearly manifests 
itself in its quantum features; the high-field magnetization shows quantum plateau at 
one third of the saturation magnetization, the susceptibility turned out to have two 
round peaks at relatively low temperatures and ESR data have proven the gap-less 
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excitation spectrum. Theoretical interest focused on the diamond chain structure also 
enhances its another experimental realizations provided by the polymeric compounds 
such as Cu 2 OS0 4 |H2| and M 3 (OH) 2 (M = Ni, Co, Mn) ^-[SH!. 

With this background, we shall investigate in the present article a simplified version 
of the frustrated Ising-Heisenberg diamond chain, which can exactly be solved by 
applying an accurate map based on the generalized decoration-iteration transformation 
|57| It is noteworthy that this relatively simple and straightforward analytical 

approach has recently been adapted to study an appearance of quantized plateaux in 
the magnetization process of the trimerized [5^1 ED] and tetramerized |!JT) Ising- 
Heisenberg linear chains. As we shall show hereafter, the same strategy can also be 
used to explore the frustrated Ising-Heisenberg diamond chain, which represents another 
particular realization of models tractable within the generalized algebraic maps. The 
main goal of present work is to exploit a grand advance of the method used for obtaining 
exact results for all possible ground states as well as all basic thermodynamic quantities. 

The rest of this paper is organized as follows. The model Hamiltonian and major 
features of the mapping method are presented in Section 2. Section 3 is divided 
into two sub-sections. In the former one, we provide rigorous results for the spin-| 
Ising-Heisenberg diamond chain, while the latter one comprises magnetic data of the 
mixed spin-^ and spin-1 Ising-Heisenberg diamond chain. Altogether, exact results for 
ground-state phase diagrams, magnetization, entropy, susceptibility and specific heat 
are derived and particularly discussed for both investigated diamond chains. A cooling 
rate of adiabatic demagnetization is also explored in connection with the enhanced 
magnetocaloric effect. Finally, some concluding remarks are drawn in Section 4. 

2. Model and its exact solution 

Let us begin by considering two kinds of spins regularly distributed on the ID lattice 
composed of the diamond-shaped units as diagrammatically depicted in Fig. 1. To 
ensure an exact tractability of this spin system, we shall further suppose that each 
diamond-shaped plaquette consists of two dumbbell Heisenberg spins (S), which are 
placed in between two Ising spins (/i) residing corner- sharing positions on the diamond 
motifs. For further convenience, we shall write the total Hamiltonian as a sum over 
plaquette Hamiltonians, i.e. fi = J2k=i^-k, where each plaquette Hamiltonian 
involves all interaction terms associated with one diamond-shaped unit (see Fig. 1) 

- FH^Ifc-i + S% k ) - Hi{p,l k _ 2 + fi z 3k+1 )/2. (1) 

Above, fx% and S k (a = x,y,z) denote spatial components of the spin-^ and spin- 
s' operators, the parameter Jh stands for the anisotropic XXZ interaction between the 
nearest-neighbouring Heisenberg (dumbbell) spins and the parameter J\ accounts for the 
Ising-type coupling between the Heisenberg spins and their nearest Ising neighbours. 
The parameter A allows to control the Heisenberg exchange interaction Jh between 
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the easy-axis (A < 1) and easy-plane (A > 1) regime and finally, the last two terms 
incorporate coupling of the Ising and Heisenberg spins to an external longitudinal 
magnetic field Hi and Hn, respectively. 

The crucial point of our calculation represents calculation of the partition function. 
By making use of the commutation rule between different plaquette Hamiltonians, i.e. 
by exploiting [H k , Hi] = valid for each k ^ I, the partition function of Ising-Heisenberg 
diamond chain can be partially factorized into the following product 

N 

Z = £ Y[Ti k ex P (-PH k ), (2) 

{(ii} k=l 

where f3 = l/{ksT), k-Q is being Boltzmann's constant, T the absolute temperature, the 
symbol means a summation over all available configurations of Ising spins {fii} 

and Tr^ stands for a trace over spin degrees of freedom of two dumbbell Heisenberg 
spins belonging to kth diamond unit. After performing this partial trace, the structure 
of relation (J2J) immediately implies a possibility of applying the generalized decoration- 
iteration mapping transformation |57| 158] 

Tr fe exp(-/%) = exp[/3F I ( / i| fe _ 2 + ^ fe+1 )/2]G[/3J I ( / i^_ 2 + ^ fc+1 )/2] 

)/2], (3) 

which converts the Ising-Heisenberg diamond chain into the uniform spin-^ Ising linear 
chain with an effective nearest-neighbour coupling R and an effective external field H . 
The expression G(x) in the first line of Eq. (J3J) depends on the spin of Heisenberg 
atoms and its explicit form is given in Appendix for two particular spin cases S = | 
and 1. It is noteworthy that a general validity of the mapping relation (JHJ) necessitates a 
self-consistency condition to be satisfied, what means, that it must hold independently 
of spin states of both Ising spins ^zk-2 and ^zk+i- It can readily be proved that a 
substitution of four possible configurations of the Ising spins ^k-2 and ii^k+i into 
the formula © indeed gives just three independent equations, which unambiguously 
determine the unknown mapping parameters A, R and Ho 

A = {G l G 2 G\) ll \ pR=]nfc*j^J t (3H = (3H l + \n(^j. (4) 

Here, we have defined the functions G\ = G(f3Ji), G2 = G(—/3Ji) and G3 = G(0) in 
order to write the transformation parameters A, R and H in more abbreviated and 
general form. Now, a direct substitution of the transformation Q into the expression 
(J2J) yields the equality 

Z(J3, Jj, J H , A, H h H B ) = A N Z ((3, R, H ), (5) 

which establishes an exact mapping relationship between the partition function Z of 
the Ising-Heisenberg diamond chain and, respectively, the partition function Zq of the 
uniform spin-^ Ising linear chain defined by means of the nearest-neighbour coupling R 
and the effective field H . Notice that the exact solution for the partition function of the 
latter system is well-known (see for instance Ref. and hence, the relation (0) can 
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readily be utilized for calculation of some important quantities (magnetization, entropy, 
specific heat, susceptibility) by making use of the standard thermodynamical-statistical 
relations. It should be mentioned, however, that there also exists another alternative 
approach that is of particular importance if some relevant physical quantity cannot be 
obtained within this procedure. Actually, the problem connected with the calculation of 
correlation functions and/or quadrupolar moment can be simply avoided by employing 
the following exact spin identities 

))o, (6) 

If ( ~z oa Q7 r, z W — l ^ Ik f 2 ^lk-2-> ffsfc-l? ^Iki frlk+l) ex P(~/^4) \ 

\ Tr k exp(-fjH k ) / 

In above, the function fx depends exclusively on the Ising spin variables {/ij}, while 
the function f 2 may depend on any spin variable belonging to fcth diamond plaquette. 
The superscripts a, 7 G (x, y, z) label spatial components of the appropriate spin 
operators and finally, the symbols (...) and (...)o stand for the standard canonical 
average performed over the ensemble defined on the Ising-Heisenberg diamond chain 
and its corresponding Ising chain, respectively. By applying the exact spin identities (JBJ) 
and (|7ji. one easily attains rigorous results for the sub-lattice magnetization m? and m| 
reduced per one Ising and Heisenberg spin, respectively, the pair correlation functions 
c hh5 c hh5 c fh5 as we U as 5 the quadrupolar moment 

m i = (PZk-2) = (^-2)0 = mo, (8) 

C U = (^-2^+1) = (/ t 3fe~2/ t 3A : +l)o = £0, (9) 

ml = (4\) = Ki/4 + m Li + e M 1: (10) 

4h = (Sk-iSk) = K 2 /A + m L 2 + e M 2 , (11) 

ell = (S x 3k -iS x 3k ) = K,/A + m L 3 + e M 3 , (12) 

4* = </&-2$a*-i> = ^1/8 + ™o(^i + MO/4 + e Li/2, (13) 

C = ((4Vi) 2 ) = ^4/4 + m L 4 + e M 4 . (14) 

As one can see, all aforelisted quantities can be expressed in terms of the single-site 
magnetization mo and the nearest-neighbour correlation Eq of the spin-| Ising linear 
chain given by R and H Q . Since exact analytical formulae for those quantities have been 
retrieved in the literature many times before [63J, we shall restrict ourselves for brevity 
to listing the coefficients K iy Li and M, (i = 1-3) emerging in the set of Eqs. ©-{El) 

Ki = FitfJi) + Fi{-(3J{) + 2^(0), (15) 
U =F i (pj 1 )-F i (-pj I ), (16) 
Mi = FiifHx) + Fti-pJx) - 2Fi(0). (17) 

An explicit representation of the functions Fi(x) is too cumbersome to write it down 
here and it is therefore left for Appendix. 

Finally, we just simply quote the well-known thermodynamical-statistical relations, 
which have been utilized for calculating Gibbs free energy Q, entropy S, specific heat C 
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and susceptibility x- Accurate results for these quantities have been obtained with the 
help of the precise mapping relation (jSJ) and using 

g = -k B T \nZ = g - Nk B T In A, (18) 

s = -{dr)H> c = - t (w^h' x = -{dip)T> (19) 

where (?o is referred to as the Gibbs free energy of corresponding spin-| Ising chain [63J. 
3. Results and Discussion 

Before proceeding to a discussion of the most interesting results it is worthy to mention 
that the results derived in the preceding section hold regardless whether ferromagnetic or 
antiferromagnetic interactions are assumed. As we are mainly interested in examination 
of the spin frustration effect, in what follows we restrict both exchange parameters 
Jh and J\ to positive values in order to match the situation in the frustrated 
antiferromagnetic diamond chain. Furthermore, it is convenient to reduce the number 
of free parameters by assuming equal g-factors for the Ising and Heisenberg spins, i.e. 
by imposing the same parameter representing the effect of external field H\ = H# = H. 
To simplify further discussion, we shall also introduce a set of reduced parameters 
t = k B T/Ji, h = H/Ji and a = Jr/Ji as describing dimensionless temperature, external 
field and the strength of frustration, respectively. 

3.1. Spin-^ Ising-Heisenberg diamond chain 

First, let us take a closer look at the ground state of the spin-| Ising-Heisenberg diamond 
chain. For illustration, two ground-state phase diagrams are displayed in Fig. 2, the 
one in the A — a space for the system without external magnetic field (Fig. 2a) and 
the other one in the a — h space for the system placed in the non-zero external field 
under the assumption A = 1.0 (Fig. 2b). Both the figures clearly demonstrate that a 
competition between the interaction parameters a, A and h gives rise to three possible 
ground states. Apart from the usual ferrimagnetic phase (FRI) and the frustrated phase 
(FRU) found both in a presence as well as absence of the external field, the system ends 
up in the saturated paramagnetic phase (SPP) once the external field is being above 
its saturation value. Spin order within FRI, FRU and SPP can be distinguished from 
one another according to their attributes (physical quantities included in the set of Eqs. 
(jKj)- (fHJ) ). as well as, through their wave functions 

N N 

|Fitf)=n|->3fc-2lII + >+>3*-l,3fc, 
k=l k=l 

ml = _0.5, m z h = 0.5, cg[ = 0.25, = 0, cf h 2 = -0.25; (20) 

N N 1 

iF^> =n i ±) 3 M n^(i+,-M-,+>) 3l _ 1 , 3i , 

k=l k=l V * 

m[ = 0, ml = 0, cf h = -0.25, = -0.25, c£ = 0; (21) 



Geometric frustration in the Ising-Heisenberg diamond chains 



7 



N N 

\spp) = ni+> 3fc - 2 ni+>+> 

k=l k=l 

ml = 0.5, ml = 0.5, = 0.25, = 0.25, c££ = 0, cf h z = 0.25. (22) 

The first product in the aforelisted eigenfunctions is carried out over all Ising spins, 
the second one runs over all pairs of the Heisenberg dumbbell spins and |±) denotes 
standard ket vector assigned to zth projection of the Ising (/i z = ±|) and Heisenberg 
(S z = ±|) spins. It is quite evident from the set of Eqs. (ffilj) that FRI displays the 
classical ferrimagnetic spin arrangement usually observed in the pure Ising systems, 
actually, all the results clearly indicate antiparallel alignment between the nearest- 
neighbouring Ising and Heisenberg spins. It should be emphasized, however, that the 
classical ferrimagnetic order originating from the antiferromagnetic Ising interaction J\ 
can be destroyed through the competing Heisenberg interaction Jh(A) that brings a 
frustration into play. As a matter of fact, the spin order dramatically changes when the 
frustration parameter a exceeds the boundary value 

(FRI\FRU) : a h = -^-p (23) 

In such case, all Heisenberg spin pairs create singlet dimers and on account of this 
singlet pairing, all Ising spins become completely free to flip. In other words, the Ising 
spins are thoroughly uncorrelated in FRU on account of the frustration arising from the 
singlet pairing between the Heisenberg spins as also suggested by (J2lf . Owing to this 
fact, it might be concluded that the diamond spin chain splits into a set of independent 
monomers (Ising spins) and dimers (Heisenberg spin pairs) whenever the frustration 
parameter is stronger than its boundary value at,. Thus, FRU is virtually being 
macroscopically degenerate monomer-dimer state with the residual entropy S ieB /3N = 
ln(2) 1 / 3 proportional to the total number of frustrated Ising spins. For completeness, 
it should be also noticed that sufficiently strong external field stabilizes the standard 
SPP regardless whether FRI or FRU constitutes the zero-field ground state (see Fig. 
2b). As could be expected, all Ising as well as Heisenberg spins tend to align into the 
external-field direction above the saturation field, which represents lower bound for an 
occurrence of SPP. 

Next, we turn our attention to the magnetization process at zero as well as non-zero 
temperatures. For this purpose, two typical magnetization vs. field dependences are 
plotted in Fig. 3 for several dimensionless temperatures. It can be readily understood 
by comparing the displayed magnetization curves with the phase diagram shown in 
Fig. 2b that two different zero-temperature limits obviously reflect both possible types 
of field-induced transitions FRI-SPP (Fig. 3a) and FRU-SPP (Fig. 3b). Since the 
ground state is being formed in the former (latter) case by FRI (FRU), the zero- 
temperature magnetization curve depicted in Fig. 3a (Fig. 3b) starts from non- 
zero (zero) magnetization in the limit of vanishing external field. Contrary to this, 
both magnetization curves always start from zero magnetization (disordered state) at 
any finite temperature according to the one- dimensional character of the investigated 
spin system. It should be emphasized, moreover, that the magnetization jumps to be 
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observed in the magnetization curves strictly at t = are gradually smeared out when 
temperature is raised from zero. In addition, the higher the temperature, the smaller 
the width of magnetization plateaux (horizontal regions in the magnetization vs. field 
dependence), which entirely disappear from the magnetization curves above a certain 
temperature. Finally, it is quite interesting to mention that the identified magnetization 
plateaux at one third of the saturation magnetization satisfy the Oshikawa-Yamanaka- 
Affleck rule [66J proposed for the formation of quantized plateaux. 

Now, let us step forward to a discussion concerned with the thermal dependence of 
the zero-field susceptibility times temperature (xt) data as displayed in Fig. 4. If the 
frustration parameter a is selected so that FRI constitutes the ground state (Fig. 4a), 
then, xt data exhibit a round minimum prior to exponential divergence appearing on 
temperature decrease. As it can be clearly seen from Fig. 4a, the stronger the frustration 
parameter a, the deeper the notable minimum whose position is simultaneously shifted 
towards lower temperatures. It should be mentioned, moreover, that an appearance of 
the round minimum in the t— xt dependence is a typical feature of quantum ferrimagnets, 
since xt product monotonously decreases with temperature for ferromagnets, while it 
monotonously increases for antiferromagnets [67J. Accordingly, a location of the round 
minimum can be regarded as a point that determines ferromagnetic-to-antiferromagnetic 
crossover in view of thermal exitations. On the other hand, if the frustration parameter 
a drives the system into the disordered FRU ground state (Fig. 4b), the xt product then 
exhibits a round minimum before it tends to the constant value 1/12 by approaching zero 
temperature. Notice that this zero-temperature value can be explained in compliance 
with the Curie law of the frustrated Ising spins, which effectively form isolated spins 
(monomers) in FRU. 

Another quantity, which is important for overall understanding of thermodynamics, 
is being the specific heat. Temperature variations of the zero- field specific heat are 
depicted in Figs. 5a-b for A = 1.0 and several values of the frustration parameter a. 
According to these plots, there still emerges at least one round maximum, which can 
be thought of as the usual Schottky-type maximum irrespective whether FRI or FRU 
constitutes the ground state. If the frustration parameter is selected sufficiently close to 
FRI-FRU phase boundary given by «b, however, there also appears an additional second 
maximum located in the low-temperature part of the specific heat. Apparently, the low- 
temperature peak becomes the more pronounced, the closer is a selected to a^. When 
the frustration parameter a drives the system into the disordered FRU ground state, 
then, the striking second maximum gradually disappears upon further strenghtening of 
a, since the low-temperature peak shifts towards higher temperatures until it entirely 
merges with the high-temperature Schottky-type maximum (see the curves labeled as 
a = 1.1, 1.25 and 1.5). These observations would suggest that the double-peak structure 
in the specific heat curves originates from thermal excitations between the ground-state 
spin configuration and the ones close enough in energy to the ground state. 

The situation becomes even more intriguing on applying the external magnetic field. 
As one can see from Fig. 5c, the rising external field generally causes a gradual increase 
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in the height of the low-temperature peak and moves it towards the higher temperatures. 
Similar trend is also seen in a change of the size and position of the high-temperature 
maximum, albeit this Schottky-like peak moves towards higher temperature more slowly 
than the low-temperature one. As a result, both maxima coalesce at a certain value 
of the external field and above this value, the specific heat exhibits just single rounded 
maximum (see for instance the curve h = 0.5). Apart from these rather trivial findings, 
a remarkable triple-peak specific heat curves can also be detected when small but non- 
zero external field is applied to the system driven by the frustration into the disordered 
FRU state (the case h = 0.05 in Fig. 5c). Besides two aforedescribed peaks, whose origin 
has been resolved earlier, there also appears an additional third peak to be located at 
lowest temperature. There are strong indications that an appearance of this additional 
sharp maximum can be explained through the field-induced splitting of the energy levels 
related to the frustrated Ising spins. In accordance with this statement, the insert of 
Fig. 5d clearly demonstrates how this peculiar third maximum gradually shifts towards 
higher temperatures with increasing the field strength until it coalesces with the second 
low-temperature peak (see the curve labeled h — 0.1). 

At last, we shall briefly discuss the adiabatic demagnetization curves studied 
in connection with the enhanced magnetocaloric effect. Some interesting results for 
adiabatic processes keeping entropy constant are presented in Fig. 6 in the form of the 
temperature vs. external field dependence. Two depicted sets of demagnetization curves 
reflect two available adiabatic scenarios related to SPP^FRI (Fig. 6a) and SPP^FRU 
(Fig. 6b) transitions. Apparently, the maximal cooling rate emerges in the vicinity of 
critical fields and zero field, where zero temperature is in principle reached whenever 
the entropy is set equal to or less than its residual value S res . It should be pointed 
out that a relatively fast heating, which occurs when the external field is lowered from 
its critical value, prevents a practical use of the whole demagnetization curves for a 
cooling purpose. In addition, the cooling effect becomes of technological relevance only 
if a cooling rate exceeds the one of paramagnetic salts. From this point of view, the 
enhanced magnetocaloric effect is observed only if the frustration drives the system into 
disordered FRU ground state (Fig. 6b). Even under this condition, the cooling rate 
of paramagnetic materials is exceeded only if the entropy is chosen close enough to its 
residual value S res and the external field is below h ~ 0.05. This limitation would imply 
that temperatures in the sub-Kelvin range are in principle accessible provided that the 
frustrated diamond chain compound with exchange constants of the order of few tents 
Kelvin, such as azurite, is used as refrigerant. 

3.2. Mixed spin-j and spin-1 Ising-Heisenberg diamond chain 

In this part, we shall turn our attention to the mixed spin-| and spin-1 Ising-Heisenberg 
diamond chain with the aim to clarify the impact of integer-valued Heisenberg spins 
on the magnetic behaviour of the frustrated diamond chain. We start our discussion 
repeatedly with the ground-state analysis. The phase diagram constructed in an absence 
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of the external field (Fig. 7a) implies existence of three possible ground states. Besides 
the semi-classically ordered ferrimagnetic phase (FRI), there also appear the quantum 
ferrimagnetic phase (QFI) and the frustrated phase (FRU). FRI, QFI and FRU can be 
characterized by means of 

N N 

\fri) = n i-> 3fc -2 n im>3*-i,3*, 

k=l k=l 

m\ = -0.5, ml = 1, qH = 1, eft = 1, eg = 0, = -0.5; (24) 



\QFI) 



N N I 

ni->s*- 2 n-^(ii.o>-io.i 

k=l k=l 



N N 

-±=(\1.0) - 10. 1^ 

3fc-l,3fc 

ml = -0.5, ml = 0.5, q** = 0.5, cf h = 0, = -0.5, cf h z = -0.25; (25) 
— ^k/tf + lfll.-l^ + I - 1.1)) - y/2^/S-l\Q.Q) 

3fc-l,3fc 



iV JV i 

|*W> = [] |±>3*-2 II 77! fv^+T(|l, -1) + I - 1, 1)) - vW^IlO, 0) 
fc=i fc=i 2 V5 L v 7 

mf = 0, m£ = 0, g£ = -c^ h = (1 + r 1 )^, c£ = -2M" 1 , eff = 0; (26) 



where 5 = y/1 + 8A 2 , the first product in the aforelisted eigenfunctions is taken over all 
Ising spins (|±) stands for fi z = ±|) and the second product runs over all pairs of the 
Heisenberg dumbbell spins (| ± 1, 0) is assigned to S z = ±1,0). Analytic expressions for 
the phase boundaries depicted in Fig. 6a read 

(QFI\FRI) : a hl = ^- [ ; (FRU\QFI) : a h2 = (27) 

The most significant difference between the two investigated diamond chains apparently 
rests in a presence of QFI located in between FRI and FRU. This observation would 
suggest that the geometric frustration initially favours uprise of QFI before it finally 
energetically stabilizes FRU. Accordingly, it might be concluded that there may not 
occur in the mixed spin-| and spin-1 diamond chain a direct frustration-induced 
transition between FRI and FRU except the one observable in the Ising limit (A = 0). 

Next, the ground-state phase diagram reflecting the effect of external field is shown 
in Fig. 7b. This phase diagram suggests that in response to the applied external 
field there also may arise the quantum ferromagnetic phase (QFO) and the saturated 
paramagnetic phase (SPP) besides the aforementioned FRI, QFI and FRU phases 

N N 



iV iV 1 

i^o) = ni+).u-2ii^(ii>o)-io,i 

fe=i fe=i v * 



3fc-l,3fc 



ml = 0.5, ml = 0.5, q& = 0.5, cf h = 0, c££ = -0.5, cf h 2 = 0.25; (28) 

N N 

3fc-l,3fcj 

fc=l fc=l 

ml = 0.5, ml = 1, q£ = 1, tfi = 1, = 0, cf h 2 = 0.5. (29) 

Before proceeding further, let us make few comments on all these possible ground states. 
Since FRI and SPP are commonly observed also in the semi-classical Ising spin systems, 
we should therefore concentrate on QFI, QFO and FRU, in which quantum entanglement 
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of the Heisenberg spin pairs indicates the quantum nature of these phases. It can easily 
be understood from Eqs. and (J28j) that QFI and QFO are quite similar one to 

each other. As a matter of fact, all pairs of Heisenberg spins reside in both phases 
the same eigenstate and thus, the only difference between them consists in the spin 
arrangement of their Ising counterparts. In QFO, which is stable at stronger fields, the 
Ising spins tend to align towards the external field, whilst they are oriented in opposite 
to the external field in QFI, which is stable at relatively weaker fields. Notice that the 
quantum entanglement between the pairs of Heisenberg spins that occurs in QFI and 
QFO can also be understood within the valence-bond-solid picture [HUES]- If spin-1 
sites are decomposed into two spin-| variables, then, one of the decomposed spins at each 
spin-1 site forms the singlet-dimer with its nearest-neighbouring spin-1 site, while the 
other one is polarized by the external field. As a result of this incomplete pairing, each 
spin-1 site effectively acts in QFI and QFO as it would be the spin-| variable. Further, 
it should be also remarked that all Ising spins are completely free to flip (frustrated) 
in FRU on behalf of a preferred antiferromagnetic alignment between each pair of the 
Heisenberg spins. Owing to this fact, FRU can be viewed as a state characterized by 
a complete randomization of the Ising spins, what consequently leads to a macroscopic 
degeneracy of FRU resembled in its residual entropy S res /3N = ln(2) 1 / 3 . 

To clarify the magnetization scenario available for the mixed spin-| and spin-1 
diamond chain, we depict in Fig. 8 all possible types of magnetization curves. In 
agreement with the phase diagram shown in Fig. 7b, there are in total four different 
types of magnetization curves reflecting the transitions FRI-SPP (a), QFI-FRI-SPP (b), 
FRU-FRI-SPP (c), and FRU-QFO-SPP (d). It should be stressed, nevertheless, that 
the system undergoes true transitions between those phases merely at zero temperature, 
where indeed emerge stepwise magnetization curves with abrupt change(s) of the 
magnetization at critical field(s). However, it is worthy to notice that there are no real 
magnetization jumps at any finite temperature and the sharp stepwise magnetization 
curves to be observed at zero temperature are gradually smeared out as temperature 
increases. It can be clearly seen from Fig. 8 that an increase in temperature actually 
shrinks also the width of plateaux until the plateau states completely disappear from 
the magnetization curves above a certain temperature. The most notable magnetization 
curves are those with the zero-field ground state corresponding to FRU as shown in Figs. 
8c-d. According to these plots, it can easily be realized that all frustrated Ising spins 
tend to align to the external-field direction for any finite but non-zero external field 
provided that temperature is set to zero. At non-zero temperatures, on the contrary, 
the magnetization rises much more steadily in the vicinity of zero field in comparison 
with the magnetization curves having the long-range-ordered FRI (Fig. 8a) and QFI 
(Fig. 8b) phases in the ground state. Finally, it is also worthwhile to remark that the 
observed magnetization plateaux at | (QFI and FRU) and | (QFO) of the saturation 
magnetization satisfy the Oshikawa-Yamanaka- Affleck rule j66j . 

Now, let us investigate in particular an influence of the spin frustration on thermal 
variations of the zero-field susceptibility times temperature (xt) data. The temperature 
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dependence of xt product is displayed in Fig. 9a for several values of the frustration 
parameter a that determine the ground state to be either of FRI or QFI type. 
Interestingly, xt data then exhibit a round minimum upon cooling, which is followed 
by an exponentially steep increase that appears by approaching zero temperature. It 
is quite obvious from Fig. 9a that the stronger the parameter of frustration a, the 
deeper the notable minimum whose position is simultaneously shifted towards lower 
temperatures. It is worthy to remember, moreover, that the minimum in the t — xt 
plot is being a typical feature of quantum ferrimagnets |HZ| and also in our case, the 
minimum becomes especially marked by selecting a G (|, 1) when QFI constitutes the 
ground state. On the other hand, the temperature dependence of xt data are depicted 
in Fig. 9b for the case when the frustration parameter a > 1 drives the system into the 
disordered FRU ground state. Under these circumstances, the susceptibility diverges 
as t' 1 at low temperatures and xt product tends to constant value 1/12 nearby zero 
temperature. This value can be interpreted in terms of the Curie law of the frustrated 
Ising spins, since the spin arrangement that appears in FRU can be viewed as being 
composed of an independent set of the antiferromagnetic Heisenberg dimers and isolated 
Ising monomers, whereas the former ones do not contribute to xt product in the zero- 
temperature limit. 

To gain an insight into overall thermodynamics, let us turn our attention to a 
thermal behaviour of the specific heat. For this purpose, the zero-field specific heat is 
plotted against temperature in Figs. lOa-b for several values of the frustration parameter 
a. It can be clearly seen from Fig. 10a that a rounded Schottky-type maximum 
observable at smaller values of a (e.g. a = 0.25) gradually changes, as the frustration 
strengthens, to a striking dependence with the shoulder superimposed on this round 
maximum (a = 0.4). It is noteworthy that the shoulder becomes the more pronounced, 
the closer is a selected to ctbi = 0.5 determining the phase boundary between FRI and 
QFI. Under further increase of the frustration the shoulder merges with the Schottky- 
type maximum what eventually gives rise to a peculiar non-rounded maximum with 
almost constant value of the specific heat over the wide temperature range (a = 0.6). 
Next, Fig. 10b shows how the specific heat recovers its double-peak structure when 
a approaches another phase boundary between QFI and FRU to appear at «b2 = 1-0 
(see for instance the curves for a = 0.75 and 0.9). In addition, it is quite evident 
from Fig. 10b that repeated strengthening of the frustration results in a suppression of 
the low-temperature peak until it finally merges with the high-temperature one. Note 
furthermore that the high-temperature peak has in general tendency to enhance in 
magnitude (both in height as well as in width) as a increases and its position is shifted 
towards higher temperatures. Altogether, it might be concluded that the remarkable 
double-peak structure of the specific heat arises just when the frustration parameter 
is close enough to a phase boundary. This result is taken to mean that the observed 
double-peak specific heat curves always originate from thermal excitations to a spin 
configuration rather close in energy to the ground state. 

Even more striking situation emerges by turning on the external field. The overall 
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trends of the external field is to increase height and width of the low-temperature peak 
and to shift it towards higher temperatures until it coalesces with the higher-temperature 
maximum. It should be remarked, nevertheless, that the most notable dependences of 
the specific heat arise from when the frustration leads to the disordered FRU ground 
state. Besides the aforedescribed general trends illustrated in Fig. 10c, there also 
appears a remarkable kind of the specific heat curve with the triple-peak structure as 
particularly drawn in Fig. lOd and its insert. Apparently, the additional third peak 
observable at very low temperatures occurs on applying the small but non-zero external 
field. It is therefore quite reasonable to conjecture that an origin of this low-temperature 
peak lies in the field-induced splitting of energy levels related to the frustration of the 
Ising spins. Actually, the stronger the external field, the greater the splitting caused 
by the external field and consequently, the position of this peak steadily shifts towards 
higher temperatures. Above certain external field, the additional third peak vanishes 
because of merging with the low-temperature peak observable also in the zero-field case. 

Finally, we shall close our discussion with an exploration of the adiabatic 
demagnetization examined in connection with the enhanced magnetocaloric effect. 
Adiabatic processes keeping the entropy constant are plotted in Fig. 11 in the form of 
the temperature vs. external field dependence. Two displayed sets of demagnetization 
curves evidently reflect adiabatic change of temperature, which accompanies the 
transitions SPP-»FRI-»QFI (Fig. 11a) and SPP^QFO^FRU (Fig. lib). When 
comparing these results with the ones formerly discussed for the spin-| diamond chain 
(Fig. 6), one easily finds some similarities between the two investigated diamond chains. 
Indeed, the most obvious drop in temperature is retrieved once again in the vicinity of 
critical fields and zero field, where zero temperature is in principle reached whenever 
the entropy is set equal to or less than its residual value S Tes . Moreover, the enhanced 
magnetocaloric effect with the cooling rate exceeding the one of paramagnetic salts 
occurs just as the disordered FRU phase constitutes the ground state and the external 
field is below h ~ 0.05. By contrast, the most obvious difference between the two 
investigated spin systems consists in the greater diversity of the adiabatic process of 
the latter (mixed-spin) system, which exhibits dependences with two critical fields in 
addition to the ones with one critical field only. 

4. Conclusion 

Exactly solvable frustrated spin systems are currently much sought after in the field of 
condensed matter physics, since they can serve as useful model systems for in-depth 
understanding of the effect of geometric frustration still not fully elucidated, yet. In 
the present work, we have provided rigorous results for one notable example of such 
system, the mixed spin-^ and spin-S* Ising-Heisenberg diamond chain tractable within 
the generalized decoration-iteration map. It is worthy to mention that this rather simple 
model system has primarily been developed to predict and to understand the behaviour 
of insulating magnetic materials, in which the Heisenberg dimers interact with the Ising 
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monomers in such a way that they form the diamond chain. Notice that the coordination 
polymers consisting of the pairs of interacting transition-metal elements (Heisenberg 
dimers) coupled to the rare-earth elements (Ising monomers) represent perspective 
experimental realizations of the proposed system. Although we are not aware of any 
real coordination compound, which would meet this requirement, the recent progress in 
the design and controlled synthesis of the molecular-based magnetic materials supports 
our hope that it would be possible to prepare such polymeric chains in the near future. 

It is worthwhile to remark that a special emphasis is in our study laid on the 
investigation of geometric frustration generated by the competition between Heisenberg- 
and Ising-type interactions. Our results clearly demonstrate that an interplay between 
the frustration and quantum fluctuations, both arising from the Heisenberg exchange 
interaction, is being at the origin of interesting behaviour not commonly observed in 
the semi-classical Ising spin systems. Indeed, this interplay gives rise to several peculiar 
ground states with entangled states of Heisenberg spins and quantum effects turned 
out to play a substantial role in determining their finite-temperature properties, as 
well. Among other matters, we have found rigorous evidence for appearance of the 
quantized plateux in the magnetization curves, the round minimum in temperature 
dependence of the susceptibility times temperature data, the double-peak zero-field 
specific heat curves, the enhanced magnetocaloric effect and so on. In our opinion, the 
most interesting finding to emerge from our study is a direct evidence of the triple-peak 
specific heat curve that appears when applying small external field to the system driven 
by the frustration into the disordered state. To the best of our knowledge, the discovery 
and possible explanation of the triple-peak structure in the specific heat curve has not 
been reported for any frustrated system hitherto. 

Last but not least, it should be remarked that the relative ease of generalized 
mapping method used here implies a possibility of further extensions. Actually, this 
approach can straightforwardly be extended to account also for the single-ion anisotropy 
effect, the biquadratic exchange interaction, the next-nearest-neighbour interaction, the 
multispin interactions, etc. It is also noteworthy that the applied procedure is not 
constrained neither by the lattice topology and thus, it can be utilized for investigating 
the effect of geometric frustration on the planar Ising-Heisenberg lattices composed of 
the diamond-shaped units [BEJEII]- In this direction continuous our next work. 

5. Appendix 

Before we list explicit expressions for the function G{x) and F^x) (i = 1-4), which enter 
into the transformation formulae ©-© and the set of Eqs. (fl5 |) - fTTj) . respectively, let 
us establish the notation j'h = /3Jn, hn = PH^ and 5 = vl + 8A 2 . 

(a) Expressions for the spin case S = |: 





G(x) = 2exp(-j H /4) cosh(x - h H ) + 2exp(j H /4) cosh(j H A/2), 
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F x (x) = - exp(-j H /4) sinh(a; - h H )/G(x), 

F 2 (x) = [expO'n/4) cosh(j H A/2) - exp(-j H /4) cosh(:r - h K )]/[2G(x)}, 
F 3 (x) = -exp(jH/4)smh(j H A/2)/[2G'(x)]. 

(b) Expressions for the spin case 5 = 1: 

G(x) = 2exp(j H /2) cosh(j H 5/2) + 2exp(-j H ) cosh[2(a: - h K )] 
+ exp(j H ) + 4cosh(:r - h H ) cosh(j H A), 

= - {2exp(-j H ) sinh[2(a; - h u )} + 2sinh(x - h H ) cosh(j H A)}/G(x), 
F 2 (x) = {2exp(-j H ) cosh[2(:r - h u )} - exp(j H ) - SC - S}/G{x), 
F 3 (x) = - {2 cosh(x - h H ) sinh(j H A) + AAS}/G(x), 
F 4 (x) = {2exp(-j H ) cosh[2(x - h H )] + exp(j H ) 

+ 2 cosh(a; - h H ) cosh(j H A) + 5C + S}/G(x), 

where S = exp(j H /2) sinh(j H ^/2)/5 and C = exp(j H /2) cosh(jn8/2)/8. 
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Figure captions 

Fig. 1 Diagrammatic representation of the mixed-spin Ising-Heisenberg diamond chain. 

The empty (filled) circles denote the lattice positions of the Ising (Heisenberg) 

spins, the ellipse demarcates fcth diamond-shaped plaquette. 
Fig. 2 (a) Ground-state phase diagram in the A — a plane for the system in an absence 

of the external magnetic field; (b) Ground-state phase diagram in the a — h plane 

for A = 1.0. 

Fig. 3 The total magnetization reduced with respect to its saturation value versus the 
external magnetic field at various temperatures t = 0.0, 0.1, 0.3, 0.5 in ascending 
order along the direction of arrows. 

Fig. 4 Thermal dependence of the zero-field susceptibility times temperature data for 
A = 1.0, a = 0.25, 0.5, 0.75, 0.9, 1.0 (Fig. 4a) and A = 1.0, a = 1.1, 1.25, 
1.5, 2.0 (Fig. 4b) in ascending order along the direction of arrows. For clarity, the 
case a c = 1.0 corresponding to the FRI-FRU phase boundary is depicted as broken 
line. 

Fig. 5 Temperature variations of the specific heat when A = 1.0 is fixed. Figs. 5a-b 
illustrate the effect of frustration parameter a on the shape of zero-field specific 
heat and Figs. 5c-d display the effect of the applied external magnetic field when 
the frustration parameter a = 1.25 drives the system into the disordered FRU state. 

Fig. 6 Adiabatic demagnetization in the form of temperature versus external field 
dependence for A = 1.0, a = 0.5 (Fig. 6a) and A = 1.0, a = 1.5 (Fig. 6b). 
For better orientation, broken curve depicts the dependence when entropy is fixed 
at the residual value S ICS /3N = ln(2)V3 Q f FRU. 

Fig. 7 (a) Ground-state phase diagram in an absence of the external field; (b) Ground- 
state phase diagram in the a — h plane for A = 1.0. 

Fig. 8 The total magnetization reduced with respect to its saturation value as a function 
of the external field at various temperatures t = 0.0, 0.1, 0.3, 0.5 in ascending order 
along the direction of arrows. 

Fig. 9 Thermal dependence of the zero-field susceptibility times temperature data for 
A = 1.0, a = 0.25, 0.4, 0.5, 0.6, 0.75, 0.9, 1.0 (Fig. 9a) and A = 1.0, a = 
1.1, 1.25, 1.5, 1.75 (Fig. 9b) in ascending order along the direction of arrows. For 
clarity, the cases a c \ = 0.5 and a C 2 = 1.0 corresponding to the phase boundaries 
between FRI-QFI and QFI-FRU, respectively, are depicted as broken lines. 

Fig. 10 Temperature variations of the specific heat when A = 1.0 is fixed. Figs. lOa-b 
illustrate the effect of frustration parameter a on the shape of zero-field specific 
heat, whereas Figs. lOc-d display the effect of applied external field when the 
frustration parameter a = 1.25 drives the system into the disordered FRU state. 

Fig. 11 Adiabatic demagnetization in the form of temperature versus external field 
dependence when A = 1.0, a = 0.75 (Fig. 11a) and A = 1.0, a = 1.75 (Fig. 
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lib). For clarity, broken curve depicts the dependence when entropy is fixed at the 
residual value S ms /3N = ln(2) 1 / 3 of FRU. 
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